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Abstract 

The origin of the pseudospin symmetry (PSS) and its breaking mechanism are explored by 
combining the supersymmetry (SUSY) quantum mechanics, perturbation theory, and similarity 
renormalization group (SRG) method. The Schrodinger equation is taken as an example, cor- 
responding to the lowest-order approximation in transforming a Dirac equation into a diagonal 
form by using the SRG. It is shown that while the spin symmetry conserving term appears in the 
single-particle Hamiltonian H, the PSS conserving term appears naturally in its SUSY partner 
Hamiltonian H. The eigenstates of Hamiltonians H and H are exactly one-to-one identical except 
the so-called intruder states. In such way, the origin of PSS deeply hidden in H can be traced in 
its SUSY partner Hamiltonian H. The perturbative nature of PSS is demonstrated by the pertur- 
bation calculations, and the PSS breaking term can be regarded as a very small perturbation on 
the exact PSS limits. A general tendency that the pseudospin-orbit splittings become smaller with 
the increasing single-particle energies can also be interpreted in an explicit way. 
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I. INTRODUCTION 



The remarkable spin-orbit splitting for the spin doublets {n,l,j = I ± 1/2), i.e., the 
spin symmetry (SS) breaking, is one of the most important concepts for understanding the 
traditional magic numbers 2, 8, 20, 28, 50, and 82 for both protons and neutrons as well 
as 126 for neutrons in stable nuclei Since these magic numbers are not simply the 

shell closure of the harmonic oscillators, it is quite sophisticated to predict the next proton 
and neutron magic numbers , which are critical for guiding the superheavy element 
synthesis. Meanwhile, the so-called pseudospin symmetry (PSS) |5|,l6| was introduced in 1969 
to explain the near degeneracy between two single-particle states with the quantum numbers 
(ra — 1, / + 2, j = / + 3/2) and (n, l,j = l + 1/2) by defining the pseudospin doublets (n = 
n — l,l = l + = /±l/2). Based on this concept, numerous phenomena in nuclear structure 
have been successfully interpreted, including superdeformation [7|, identical bands 

a a. 



and pseudospin partner bands 
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ll|, etc. It will be quite interesting and challenging to 



understand the shell closure and the pseudospin symmetry on the same footing, in particular 
for superheavy and exotic nuclei near the limit of nucleus existence. 

During the past decades, more and more exotic nuclei have became accessible with world- 
wide new generation radioactive ion beam (RIB) facilities. It has been shown that the 



traditional magic numbers can change in the nuclei far away from the stability line |12 |. 
The splittings of both spin and pseudospin doublets play critical roles in the shell structure 
evolutions, for example, the = 28 shell closure disappears due to the quenching of the 



spin-orbit splitting for the ulf spin doublets jl3l-ll6|. the Z = 64 sub-shell c 



osure is closely 



related to the restoration of PSS for the 7r2^ and ttI/ pseudospin doubltes |17H19|. There- 
fore, it is a fundamental task to explore the origin of SS and PSS, as well as the mechanism 
of their breaking in both stable and exotic nuclei. 

Since the suggestion of PSS in atomic nuclei, there have been comprehensive efforts to 
understand its origin. Apart from the formal relabeling of quantum numbers as shown 
above, various of explicit transformations from the normal scheme to the pseudospin scheme 



have been proposed 20N22|. Based on the single-particle Hamiltonian of the oscillator shell 



model, the origin of PSS is connected with the special ratio in the strengths of the spin- 
orbit and orbit-orbit interactions |20 |. The relation between the PSS and the relativistic 



mean-field (RMF) theory 



23 



25| was first noted in Ref. 26|, where the RMF theory 
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explains approximately such special ratio in the strength of the spin-orbit and orbit-orbit 
interactions. 

As a substantial progress, the PSS was shown to be a symmetry of the Dirac Hamiltonian, 
where the pseudo-orbital angular momentum / is nothing but the orbital angular momentum 
of the lower component of the Dirac spinor j^^]. In addition, the equality in magnitude but 
difference in sign of the scalar potential S{r) and vector potential V(r) was suggested as 
the exact PSS limit by reducing the Dirac equation to a Schrodinger-like equation j2^. A 



more general condition d{S + V)/dr = was proposec 
satisfied in exotic nuclei with highly diffuse potentials 



281 1. which can be approximately 
291. However, since there exist no 



bound nuclei within such PSS limit, the non-perturbative nature of PSS in realistic nuclei 



has been suggested |30| 
dynamical symmetry 



3l| . which was also related to the consideration of the PSS as a 



32|. In this sense, an explicit and quantitative connection between 
the ideal PSS limits and the realistic nuclei is still missing. 

After the PSS was revealed as a relativistic symmetry, numerous efforts have been ded- 
icated to tracing the relativistic origin of PSS and its breaking mechanism in a quanti- 
tative way. These investigations include the one-dimensional Woods-Saxon potential 33| . 



and Yukawa 

oscillator [78|, Hartmann 
some formal studies 183 



Poschl- Teller |6l 464j. diatomic molecular 65|, Manning- Rosen 66l- l68l | , M ie-type 



the spherical Woods-Saxon 



monic oscillator 49| . Hulthen [50|, |5l| , Mor se l52 



34|-|38 |. Coulomb 39l. l40l| . harmonic oscillator 
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41-48 



anhar- 



Rosen-Morse 



57, 



Eckart 



59| 

69 



60|, 
7q|, 



7l| potentials, as well as the deformed harmonic oscillator 72l-l77l| . anharmonic 



79 



80| . Kratzer [8l|, and Makarov 82[ poten tials , together with 



103 



-1871]. Self-consistently, the PSS in spherical |88l - ll02| and deformed 



109( 1 nuclei have been investigated in the RMF and relativistic Hartree-Fock 



110 



theories. The PSS investigations were also extended to the single-particle resonances 



112| 



113- 



llq as well as t he single-particle states in the Dirac sea, i.e ., the SS in the anti-nucleon 



spectra 117l - ll21| and the A spectra in hyper- nucle i 122Hl24i |. The relevances of the PSS 



in nuclear magnetic moment s and tra nsitions 125 
nucleon-n ucleon scatterings 127 - 



1261 as well as in nucleon-nucleus and 



1311. 



130[ | were discussed. The readers are referred to Refs. 



132| for reviews 



Recently, in Refs. 



133 



134j . the perturbation theory was used to investigate the symme- 



tries of the Dirac Hamiltonian and their breaking in realistic nuclei. This provides a clear 
way for investigating the perturbative nature of PSS. It is found that the energy splitting of 
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the pseudospin doublets can be regarded as a result of perturbation on the Hamiltonian with 
relativistic harmonic oscillator potentials, where the pseudospin doublets are degenerate. 
The supersymmetric (SUSY ) quantu m mechanics was also used to investigate the symme- 



tries of the Dirac Hamiltonian 



135 



5l- |l42| ■ In particular, by employing both exact and broken 
SUSY, the phenomenon that all states with / > have their own pseudospin partners except 
the so-called intruder states can be interpreted naturally within a unified scheme. A PSS 



breaking potential without singularity can also be obtained with the SUSY technique |136 |. 
In contrast, singularities appear in the reduction of the Dirac equation to a Schrodinger- 
like equation for the lower component of the Dirac spinor. However, by reducing the Dirac 
equation to a Schrodinger-like equation for the upper component, the effective Hamilto- 



nian shown in Ref. 
calculations. 



136| is not Hermitian, which prevents the quantitative perturbation 



A very recent work 1431 ] filled the gap between the perturbation calculations and the 
SUSY descriptions by using the similarity renormalization group (SRG) technique to trans- 
form the Dirac Hamiltonian into a diagonal form. The effective Hamiltonian expanding in a 
series of 1/M in the Schrodinger-like equation thus obtained is Hermitian. This makes the 
perturbation calculations possible. Therefore, we deem it promising to understand the PSS 
and its breaking mechanism in a fully quantitative way by combining the SRG technique, 
SUSY quantum mechanics, and perturbation theory. 

By using the SRG technique, a Dirac equation can be transformed into a diagonal form in 
a series of 1/M, and its lowest-order approximation corresponds to a Schrodinger equation. 
By taking this lowest-order approximation as an example, the idea for applying the SUSY 
quantum mechanics to trace the origin of the PSS will be illustrated and the PSS breaking 
mechanism will be explored quantitatively by the perturbation theory in this paper. In 
Sec. [ni the SUSY quantum mechanics will be briefly recalled with its application for the ra- 
dial Schrodinger equation. The numerical details for solving the radial Schrodinger equation 
in coordinate space and the results for the single-particle eigenstates, the pseudospin-orbit 
splittings, the superpotentials, the PSS conserving and breaking potentials, as well as the 
perturbation corrections to the single-particle energies will be presented in Sec. IIIII Finally, 
summary and perspectives will be given in Sec. IIVI 
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II. THEORETICAL FRAMEWORK 



A. Similarity renormalization group for Dirac Hamiltonian 

Within the relativistic scheme, the Dirac Hamiltonian for nucleons reads 

HD = (y-p + f3{M + S)+V, (1) 

where a and /3 are the Dirac matrices, M the mass of nucleon, and S and V the scalar 
and vector potentials, respectively. According to the commutation and anti-commutation 
relations with respect to /3, the Dirac Hamiltonian can be separated into the diagonal e 
and off-diagonal o parts, Hd = e + o, which satisfy [e, (3] = and {o, /3} = 0. In order to 
obtain the equivalent Schrodinger-like equation for nucleons, the main task is to decouple the 
eigenvalue equations for the upper and lower components of the Dirac spinors. A possible 
way is to make the off-diagonal part of the Dirac Hamiltonian vanish with a proper unitary 
transformation. 



According to the similarity renormalization group (SRG) technique 144l . Il45l |. the Hamil- 
tonian Ho is transformed by an unitary operator U (/) as 

Hd{1) = U{1)HdU\1) (2) 

with H£,{1) = e{l) + o(/), Hu{0) = H^,, and a flow parameter /. Then, the so-called 
Hamiltonian flow equation can be obtained by taking the differential of the above equation, 
i.e., 

f^Ho{l) = [v{l),Hn{l)] (3) 



with the anti-Hermitian generator rj(l) = As pointed out in Ref. 145 |. one of 



the proper choices of r]{l) for letting off-diagonal part o(/) = as / — )■ oo reads r]{l) = 
[PM, Hd{1)]. Then, the diagonal part of the Dirac Hamiltonian e{l) at the / oo limit can 
be derived analytically in a series of 1/M, 

£(00) = Meo(oo) +ei(oo) + ; , ; , 



M M2 

+ ^[[a-p,WS + V)],a-p] + --- (4) 
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In such way, the eigenequations for the upper and lower components of the Dirac spinors 
are completely decoupled. The equivalent Schordinger-like equations for nucleons with Her- 
mitia n eff ective Hamiltonians can be obtained. The corresponding details can be found in 



Ref. 



143| 



For the nucleons in the Fermi sea, the eigenequations with the Hamiltonian in Eq. 
up to the 1/M order correspond to Schordinger equations 

1 



+ Vir) 



= Eaipa{r), (5) 



2M 

where V{r) = S{r) + V(r), and the rest mass of nucleon M is reduced in the single-particle 
energies E. By assuming the spherical symmetry, the radial equations can be cast in the 
form of 

HR,{r) = EMr) (6) 

with the single-particle Hamiltonian 

H = '^ + lA!l^±}l + V(r), (7) 

and the single-particle wave functions 

Mr) = ^KUr). (8) 

where ^j^^a spherical harmonics spinors, the single-particle eigenstates are specified 

by the set of quantum numbers a = (a, nia) = {ua, la, ja, ^a), and the good quantum number 
K = + 1 /2) for j = / ± 1 /2 is adopted. 

It can be clearly seen that H conserves the explicit spin symmetry (SS) for the spin 
doublets a and b with + Kf, = —1, which leads to the same centrifugal barrier (CB) k(k + 
l)/(2Mr^). Similarly, in order to investigate the origin of the pseudospin symmetry (PSS) 
and its breaking mechanism, it is crucial to identify the corresponding term proportional to 
/(/ + 1) = k,{k — 1), which leads to the same k(k — 1) values for the pseudospin doublets a 



and b with + = As noted in Ref. 



136| . we also consider the supersymmetric (SUSY) 



quantum mechanics is one of the most promising approaches for identifying such k(k — 1) 
structure. 

In the following, we will briefly recall the key formalism of SUSY quantum mechanics 



136 



146| . Then we will focus on the application of the SUSY quantum mechanics to the 



Schrodinger equations. 



B. Supersymmetric Quantum Mechanics 



It has been shown that every second-order Ham iltonian can be factorized in a product of 
two Hermitian conjugate first-order operators jl47 |. i.e., 



= B+B- 



(9) 



with B = [-B^]^ Its SUSY partner Hamiltonian can thus be constructed by 146 1 



H. = B-B+. 



(10) 



Hi 



Since the extended SUSY Hamiltonian 



is the square of the Hermitian operators 



all eigenvalues Es{n) of the eigenvalue equation 



H2 





are non-negative, and the two-component wave functions read 




(11) 



(12) 



(13) 



(14) 



where ipi{n) and ip2{n) are the eigenfunctions of Hi and H2, respectively. It can be easily 
seen that, for each eigenstate with Es{n) > 0, it is an eigenstate for both Hi and and 
the corresponding eigenfunctions satisfy 

B- B^ 



Es(n 



--ipi{n), ilJi{n) 



ip2{n) 



(15) 



with the normalization factors l/\/EQ{n). 



The property of SUSY can be either exact (also called unbroken) or broken 146| . On 



one hand, the SUSY is exact when the eigenvalue equation (1131) has a zero energy eigenstate 
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— — 



1 2 

exact SUSY 



1 2 

broken SUSY 



FIG. 1: (Color online) Schematic patterns of the exact and broken supersymmetries (SUSY). 

Es{0) = 0. In this case, as an usual convention, the Hamiltonian Hi has an additional 
eigenstate at zero energy that does not appear in its partner Hamiltonian H2, since 



B-iji{0) = 0, ^2(0) = 0, 



(16) 



i.e., the trivial eigenfunction of H2 identically equals to zero. On the other hand, the SUSY 
is broken when the eigenvalue equation ( fT3|) does not have any zero energy eigenstate. In this 
case, the partner Hamiltonians Hi and H2 have identical spectra. The schematic patterns 
of both cases are illustrated in Fig. [T] 



C. SUSY quantum mechanics for SchrSdinger equations 

For applying the SUSY quantum mechanics to the Schrodinger equations shown in Eq. 
first of all, one sets a couple of Hermitian conjugate first-order operators as 



b: 



V2M' 



B- 



V2M 



(17) 



whe re the as the function of r are the so-called superpotentials to be determined 1361 . 
1461]. Then, the SUSY partner Hamiltonians Hi and H2 can be constructed as 



Hi{k) 

H2(k) 



B-Bt 



1 


■ <P 


2M 




1 


'' d^ 


2M 


dr"^ 



+ QI + Q'. 



(18a) 
(18b) 



In order to explicitly identify the k (k + 1) structure shown in Eq. ([7]), the reduced superpo 
tentials q^ir) are assumed as [136 1 

K 



In such way, the Hamiltonians Hi and H2 can be further expressed as 
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(19) 



BtB- 
B-Bt 



2M 
1 

2M 



d'^ k(k + 1) n 2k , 
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r> 2 

r 



(20a) 
(20b) 



It can be seen that not only the k(k + 1) structure appears in Hi but also the k,{k, — 1) 
structure explicitly appears in the SUSY partner Hamiltonian H2. The so-called pseudo- 
centrifugal barrier (PCB) term k{k, — l)/(2Mr^) leads to the conservation of the PSS. 

In general, the Hamiltonian H in the Schrodinger equation ([6]) differs from the SUSY 
Hamiltonian Hi in Eq. (!20l) by a constant, i.e.. 



H{k) = Hi{K) + e{K) 



(21) 



where e(/t) is the so-called energy shift 1361 . Il46| . The K-dependent energy shifts can be 
determined in the following ways: 1) For the case of k < 0, it is known that the most deeply 
bound state for a given k, e.g, lsi/2, 1^3/2? etc., has no pseudospin partner. This indicates 
the exact SUSY is achieved, and requires 



e K 



Eitt- 



(22) 



2) For the case of k > 0, each single-particle state has its own pseudospin partner. This 
indicates the SUSY is broken, and thus the corresponding energy shift can be, in principle, 
any number which makes the whole set of Hi eigenstates positive. In practice, the energy 
shifts in this case are determined by assuming that the pseudospin-orbit (PSO) potentials 
vanish as r — )■ 0, which is similar to the usual surface-peak type spin-orbit (SO) potentials. 

In order to fulfill the above condition, we first analyze the asymptotic behaviors of the 
reduced superpotentials ^^(r). Combining Eqs. ([7]), (]20ap . and (!2T|) . one has 



2M 



r 



+ e K 



V{r). 



At large radius, for potential lim^^oo^('") = 0, qK{r) becomes a constant as 



lim q^{r) 



-2Me{K). 



At small radius, for any regular potential V{r), it requires ^^(0) = 0, and also 



(23) 



(24) 



(25) 



ClS Sb linear function of r. 

As the PSO potentials vanish at the original poin t, it requires that limj—^o quaif') 



liiRr^Q qf^^{r) with + ft^ = 1 for pseudospin doublets [l36|. Finally, the energy shifts are 
determined by 

e(fi:,) = 2 Vl^^ - ei^b) (26) 

for the case of > 0. 

Before ending this section, it is interesting to seek a possible exact PSS limit analytically. 
First of all, combining Eqs. ([7]), f l20bp . and ( ETjl . the SUSY partner Hamiltonian of H{k) 



reads 

//W = ff,W+eW = -^ + ^|i^ + V;W (27) 

with 

K{r) = Vir) + q'^{r)/M. (28) 

In this paper, we use a tilde to denote the operators, potentials, and wave functions belonging 
to the representation of H. Then, by definition, the exact PSS limit holds EnKa = -£'(n-i)Ki, 
with < and + = 1. This indicates 

H2{Ka) + e{Ka) = HiiK^) + e(Kb). (29) 

By combining Eqs. fl20l) and (l23l) . as well as the boundary condition ^'^(0) = 0, one can 
readily have 

(lKa{r) = q^,{r) = Mu^r (30) 

with a known constant = (e(K(j) — e{K,b))/ {Kb — i^a)- As the reduced superpotentials q^ir) 
are linear functions of r, the central potential V{r) in H has the form of 

Vno{r) = ^u;y + V{0). (31) 

The corresponding PSS limit is nothing but the well known case with harmonic oscillator 
(HO) potentials, which leads to the energy degeneracy of the whole major shell. 

III. RESULTS AND DISCUSSION 

To perform a quantitative investigation, the mass of nucleon is taken as M = 939.0 MeV, 
and the central potential V{r) in Eq. ([7]) is chosen as a Woods-Saxon form 
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FIG. 3: (Color online) Pseudospin-orbit splittings (E'j^ — Ej^)/(2l + 1) versus the average single- 
particle energies {Ej^ + Ej-^)/2 for pseudospin doublets. 

with the parameters Vq = —63.297 MeV, R = 6.278 fm, and a = 0.615 fm, which corresponds 



to the neutron potential provided in Ref. 



148 



by taking = 82 and Z = 50. This potential 



is shown as the solid line in Fig. E T he radial Schrodinger equations are solved in coordinate 
space by the shooting method |l49| within a spherical box with radius -Rbox = 20 fm and 
mesh size dr = 0.05 fm. 
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6 

r(fm) 
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FIG. 4: (Color online) Single-particle wave functions Rnij{r) of H for the 3si/2, 2^3/2, 2^5/2, and 
157/2 states. 

A. Representation of single-particle Hamiltonian H 

In Fig. O the discrete single-particle states obtained in the Woods-Saxon potential are 
shown. In order to see the K-dependence and the energy-dependence of pseudospin-orbital 
(PSO) splittings more clearly, we plot the reduced PSO splittings {Ej^ — Ej^)/{21 + 1) 
versus their average single-particle energies E.^^ = {Ej^ + Ej^)/2 in Fig. |3l where j< (j>) 
denotes the states with j = / — 1/2 (j = / + 1/2). It is seen that the amplitudes of the 
reduced PSO splittings are less than 1 MeV. Moreover, as a general tendency, the splittings 
become smaller with the increasing single-particle energies. Such energy- dependent behavior 
was also found in the self-consistent relativistic continuum Hartree-Bogoliubov (RCHB) 



calculations 
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resonance states 



29 1 ■ It was also reported that the PSO splittings can even reverse in 



113 



115| . Note that the spin-orbit (SO) splittings never reverse in realistic 



nuclei. Therefore, it is very interesting to investigate the physical mechanism for such 
energy- dependent behavior. This will be also helpful to figure out whether the pseudospin 
symmetry (PSS) is an accidental symmetry 46| or not. 

In Fig. m the corresponding single-particle radial wave functions Rnij{r) of H are shown 
by taking the 3si/2, 2^3/2, 2^5/2, and 1(77/2 states, i.e., the pseudospin doublets 2p and 1/, as 
examples. In the following discussions, these two pairs of pseudospin doublets will be often 
used for illustration. It is clear that the wave functions of the spin doublets are identical 
since there is no spin-orbit term in H. In contrast, the wave functions of the pseudospin 
doublets are very different from each other. This leads to a difficulty to analyze the origin 
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TABLE I: Contributions from the kinetic term (kin.), centrifugal barrier (CB), and central potential 
(cen.) to the single-particle energies E and the corresponding pseudospin-orbit splittings AE'pso 
for the pseudospin doublets 2p and If. All units are in MeV. 



state i?kin. EqB Eccn. E 

3sy2 28.953 0.000 -50.545 -21.591 

2d3/2 16.845 11.758 -51.746 -23.143 

A^Jpso 12.109 -11.758 1.201 1.552 

2(i5/2 16.845 11.758 -51.746 -23.143 

1^7/2 6.197 20.483 -54.188 -27.508 

A^PSo 10.648 -8.725 2.442 4.365 



of the PSS and its breaking. 

Prior to the quantitative analysis by using the perturbation theory in Ref. 



133) . the 



investigation of PSO splittings Ai^pso was usually done by decomposing the contributions 
by terms, where each contribution is calculated with the corresponding operator Oi by 



R{r)OiR{r)dr. 



(33) 



Within the representation of H shown in Eq. ([7]), the operators of the kinetic term, cen- 
trifugal barrier (CB), and central potential read —d'^/{2Mdr'^), k,{k + l)/(2Mr^), and V{r), 
respectively. In Table [H the contributions from these terms to the single-particle energies E 
as well as the corresponding PSO splittings AEp^o ^i-rc shown for the pseudospin doublets 
2p and 1/. It is not surprising that, within this representation, the contributions to AEp^o 
come from all channels, while they substantially cancel to each other in a sophisticated way. 

In previous studies, the phenomenon of such s tron g cancellations among different terms 



was usually associated with the dynamical 32 



46 



143| and even the nonperturbative 30 



31 



12l| nature of PSS. However, such connection is mystified and sometimes even misleading. 
We will demonstrate by using the perturbation theory that the nature of PSS in the present 
investigation is indeed perturbative. 



The main idea for using the Rayleigh-Schrodinger perturbation theory 150| to investigate 
the spin symmetry (SS) and PSS in single-particle Hamiltonian as well as their breaking in 
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2 4 6 8 10 12 

r(fm) 



FIG. 5: (Color online) Woods-Saxon potential in H (solid line) and harmonic oscillator potential 
m H^'-' (dashed line) as a function of r. 



atomic nuclei can be found in Ref. 133| . Following this idea, the Hamiltonian H is split as 



H = Ho + W, 



(34) 



^ 1 for m ^ k 



(35) 



where Hq conserves the exact PSS and W is identified as the corresponding symmetry 
breaking potential. The condition 

Ek — E„ 

with Wmk = (^m| l^fc) determines whether W ca n be treated as a small perturbation and 
governs the convergence of the perturbation series 150 1. 

For the present case, it has been analytically shown in Section III CI that the Hamiltonian 
with harmonic oscillator (HO) potentials is one of the exact PSS limits. Thus, one has 



HO 



2M 



rf2 k(k + 1) 



f^V + no), 



(36) 



and W^'-' is just the difference between H and H^'^. To minimize the perturbations to the 
sdg states, the coefficient uj is chosen as 1.118 x 41A~^/^ MeV, and the trivial constant V^(0) 
is taken as —73 MeV, as illustrated in Fig. |5l Although the symmetry breaking potential 
pyHO (diverges at r — >■ oo due to the parabolic behavior of H^*^, the property that the bound 
state wave functions decay exponentially at large radius leads to convergent results of the 
matrix elements Wmk- 

In the upper panel of Fig. [6l the values of \W^^ / {Em — -Efc)] for the pseudospin dou- 
blets k = 2p and 1/ are shown as functions of the energy differences Em — E^, where the 
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FIG. 6: (Color online) Panel (a): Values of | W^°/(^m -^fc)| versus the energy differences — 
for the pseudospin doublets k = 2p and 1/. Panel (b): Corresponding single-particle energies 
obtained at the exact PSS limit H^^'^, and by the first-, second-, and third-order perturbation 
calculations, as well as those obtained with H. 

unperturbed eigenstates in perturbation calculations are chosen as those of H^'~'. Since 
the spherical symmetry is assumed, only the single-particle states m and k with the same 
quan tum number k lead to non- vanishing matrix elements Wmk- Similar to the results in 



Ref. 1331 ] ■ the values of \ W^^ / {Em — -Efc)| decrease as a general tendency when the energy 
differences Em — Ek increase. From the mathematical point of view, this property provides 
natural cut-offs of the single-particle states in perturbation calculations. For example, the 
first-, second-, and third-order perturbation corrections to the single-particle energies in the 
present calculations are of 0.001 MeV accuracy when the energy cut-off for Em is taken as 
150 MeV. Furthermore, it is shown that the largest perturbation correction is roughly 0.13. 
This indicates that the criterion in Eq. fl35|) can be fulfilled. 

In the lower panel of Fig. [6l the single-particle energies obtained at the exact PSS limit 
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FIG. 7: (Color online) Reduced superpotentials q^ir) for the p and / blocks. 

H^^, and their counterparts obtained by the first-, second-, and third-order perturbation 
calculations, as well as those obtained with H are shown from left to right. Although the 
perturbation corrections do not converge very fast since the largest perturbations for all 
four states are beyond 0.1, the PSO splittings are well reproduced up to the third-order 
perturbation calculations. 

The above results have pinned down the perturbative nature of PSS in the present in- 
vestigation. In the next subsection, we will study the origin of the PSS and its breaking 
mechanism in an explicit way within the representation of the supersymmetric (SUSY) part- 
ner Hamiltonian H. 

B. Representation of SUSY partner Hamiltonian H 

In order to obtain the SUSY partner Hamiltonian H in Eq. ( l27l) . first of all, one should 
solve the first-order differential equation ( l23l) for the reduced superpotentials qKir) with the 
boundary condition ^^(O) = 0. By taking the p and / blocks as examples, the corresponding 
quir) are shown in the unit of MeV/c in Fig. [71 Since the L.H.S. of Eq. f l23|) contains 
a K-dependent term oc n/r, the reduced superpotentials qui^r) thus obtained are also k- 
dependent. In contrast, it should be emphasized that qi^{r) does not depend on the radial 
quantum number n for a given k. One will discover in the following that such n-independent 
property is essential for understanding the general pattern of Ai^pgo versus Eg,^ as shown in 
Fig. |3l In addition, it can also be examined that the reduced superpotentials q^ir) satisfy 
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FIG. 8: (Color online) Upper panels: K-dependent central potentials ^^(r) in H as a function of 
r for the p and / blocks, while the Woods-Saxon potential in H is shown for comparison. Lower 
panels: The corresponding single-particle energies obtained with H and H. 

their asymptotic behaviors at r — )■ and r — )■ oo in Eqs. (125|) and (l24l) . 

The K-dependent central potentials in H can be then calculated by Eq. f l28|) . and 

the corresponding asymptotic behaviors read 

and 

lim V;(r) = 0. (38) 

1 — >oo 

It is important that these potentials are regular and converge at both r — )■ and r — )■ oo. 

In the upper panels of Fig. [SI the central potentials V^ir) in H are shown by taking the 
p and / blocks as examples, while the Woods-Saxon potential V{r) in H is also shown for 
comparison. For all blocks, the potentials K;(r) approximately remain a Woods-Saxon shape, 
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FIG. 9: (Color online) Single-particle wave functions R^lji^) of H for the 2pi/2 and If states. 

and they are shallower than the original potential V{r). Focusing on a pair of pseudospin 
partners, the potential V^ir) with k < is higher than that of its pseudospin partner with 
K > at r < 6 fm, since in this region the superpotential qnir) with k < increases faster 
with r as shown in Fig. [71 In contrast, these two potentials reverse at r > 6 fm, since 
in this region the slope of ^^(r) with k < becomes negative whereas that with k < 
remains positive. By comparing the panels (a) and (b), it is seen that the amplitude of the 
difference between Vi^{r) for a pair of pseudospin partners increases with the difference of 
their quantum numbers — ^bl- 

After getting the central potentials K;(r), it is ready to calculate the single-particle en- 
ergies and wave functions of the SUSY partner Hamiltonians H{k). In the lower panels of 
Fig. [HI the discrete single-particle energies obtained with H are shown and compared with 
those obtained with H . It is clearly shown that the eigenstates of Hamiltonians H and 
H are exactly one-to-one identical, except the lowest eigenstates with k < in iJ, which 
are the so-called intruder states. In other words, the fact that the intruder states have no 
pseudospin partners can be interpreted as a natural result of the exact SUSY for k < and 
broken SUSY for k > 0. 

By holding the one-to-one mapping relation in the two sets of spectra, the origin of PSS, 
which is deeply hidden in the Hamiltonian H within the normal scheme, can be now traced 
by employing its SUSY partner Hamiltonian H . 

In order to shed more light on this idea, we first show the single-particle radial wave 
functions of H for the 2p and 1/ pseudospin doublets in Fig. [9] In the SUSY 
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TABLE II: Contributions from kinetic term (kin.), pseudo-centrifugal barrier (PCB), and central 
potential (cen.) to the single-particle energies E and the corresponding pseudospin-orbit splittings 
Aii^PSo for the pseudospin doublets 2p and If. All units are in MeV. 



state 


-Skin. 


EpCB 


E 

-'-'ccn. 


E 


2Pl/2 


16.602 


6.723 


-44.916 


-21.591 


2P3/2 


17.331 


6.857 


-47.332 


-23.143 


^Epso 


-0.729 


-0.134 


2.415 


1.552 


1/5/2 


5.710 


16.286 


-45.139 


-23.143 


1/7/2 


6.293 


16.591 


-50.392 


-27.508 


^-E'pso 


-0.584 


-0.305 


5.253 


4.365 



quantum mechanics 146| . the numbers of nodes in the radial wave functions R^ijir) are one 
less than those in their counterparts Rnij{r) when the SUSY is exact, while the numbers 
of nodes in R^ijir) are the same as those in their counterparts Rnijir) when the SUSY is 



broken. This indicates the well known node relation between the pseudospin doublets |84| . 



n = n — 1 for k < 0, n = n for k > 0, (39) 

is nothing but an intrinsic property of the SUSY quantum mechanics. This node relation 
can be also checked by comparing the wave functions shown in Fig. |9] and Fig. HI In fact, not 
only the numbers of nodes are equal, but also the wave functions of pseudospin doublets are 
almost identical to each other. Therefore, within this representation, the quasi-degeneracy 
of pseudospin doublets is closely related to the similarity of their wave functions, and vice 
versa. 

The same strategy is then adopted to investigate the PSO splittings AiJpso by decompos- 
ing the contributions from each term as done in Table HI but now within the representation 
of H shown in Eq. f l27|l instead. The corresponding operators include the kinetic term 
—d^/(2Mdr'^), the pseudo-centrifugal barrier (PCB) k,{k — l)/(2Mr^), and the central po- 
tential Vi^{r). The corresponding results for the pseudospin doublets 2p and If are listed in 
Table Hll It can be seen that for each pair of pseudospin doublets the energy contributions 
from the PSS conserving terms, i.e., the kinetic term and PCB, are very similar. The PSO 
splittings Ai^pso are mainly contributed by the difference in the central potentials AE^cen., 
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FIG. 10: (Color online) Pseudospin symmetry conserving potentials VpssC?^) (a) and breaking 
potentials Vpsol^') (b) for the p and / blocks. 

which is due to the shght ^-dependence of V^{r) as shown in Fig. [HI In other words, the 
sophisticated cancellations among different terms in H can be clearly understood by using 
a proper decomposition with the help of the SUSY quantum mechanics. 

In order to perform the quantitative perturbation calculations, the Hamiltonian H is split 

as 

H = H^^^ + lyPS^ (40) 

where Hq^^ and ly^^^ are the corresponding PSS conserving and breaking terms, respec- 
tively. By requiring that 14^^^^ should be proportional to k, which is similar to the case of 
the spin-orbit term in the normal scheme, one has 

K„(r) = Vpss{r) + KaVpso{r), (41a) 
K,(r) = Vpss{r) + K,Vpso{r), (41b) 
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for a pair of pseudospin doublets with = 1. In such way, the PSS conserving potentials 

^ssl'") and breaking potentials Vpso('") can be uniquely determined, and the corresponding 



HP^ and l^PSS ^g^^ 



H, 



PSS 



° 2M 



, 2 

PSS 



+ V^.ss(r-), (42a) 



H^^^^ = /s:Vpso(r). (42b) 

By taking the p and / blocks as examples, the PSS conserving potentials Vpss('") and 
breaking potentials Vpso(^) are illustrated in Fig. [TOl On one hand, it can be seen that 
the PSS conserving potentials Vpss('") remain an approximate Woods-Saxon shape, and 
they are K-dependent to a small extent. On the other hand, the PSS breaking potentials 
^so(^) show several special features. First of all, these PSS breaking potentials are regular 
functions of r, in particul ar, th ey vanish at r — )■ oo. This was also one of the main goals of 
the investigations in Ref. 136| . but here we not only achieve the goal, but also keep every 
operator Hermitian. Second, it can be seen that the amplitudes of V^so are around 1 MeV, 
which directly lead to the reduced PSO splittings AE'pso ^ 1 MeV as shown in Fig. |3l 
More importantly, different from the usual SO potentials with a surface-peak shape, the 
PSO potentials Vpso('") are negative at small radius but positive at large radius with a node 
at the surface region. This property can be traced back to the differential equation ( l23l) of 
the reduced superpotentials ^'^(r). By analyzing the corresponding asymptotic behaviors of 
q^ir) at r — 7- and r — )■ oo, one can conclude that such particular shape of Vpso('^) holds as 
long as the central potentials V{r) in the Schrodinger equations being a Woods-Saxon-like 
shape. 

The particular shape of the PSO potentials Vpso('") can explain well the variations of 
the PSO splitting with the single-particle energy. First of all, it has been emphasized above 
that Vpso(?") do not depend on the radial quantum number n. Meanwhile, the single-particle 
wave functions R{r) extend to larger distance with higher energies. Thus, the matrix element 
{R\Vy>so\R) is negative when the wave function is centralized in the inner part. As the wave 
function becomes more extended, the positive part of Vpso('^) compensates the negative value 
of the matrix element. In such way, the PSO splittings Ai^pso decrease while the radial 
quantum numbers n increase. The splittings can even reverse for the resonance states, where 
the outer part of the PSO potentials plays the major role. 

Finally, we preform the perturbation calculations based on the pseudospin symmetric 
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Hamiltonian Hq^^ with the perturbation W^^^. In Fig. HH the values of W^l^/{Em - Ek 
for the pseudospin doublets k = 2p and 1/ are shown as functions of the energy differences 
Em — Ek in the upper panel, while the single-particle energies obtained at the PSS limit 
Hq^^, and their counterparts obtained by the first-, second-, and third-order perturbation 
calculations, as well as those obtained with H are shown from left to right in the lower panel. 
It can be seen that the pseudospin doublets are exact degenerate at the PSS limit Hq^^- It 
can also be examined that the radial wave functions for each pair of pseudospin doublets 
are identical. For the present decomposition, the largest perturbation correction is less than 
0.03, which is almost one order of magnitude smaller than that with the decomposition 
H = H^'~' + This indicates that the criterion in Eq. f l35|) is satisfied quite well. As 

shown in the lower panel, not only the PSO splittings but also the energy degeneracy of the 
spin doublets are excellently reproduced by the first-order perturbation calculations. 

In such explicit and quantitative way, the PSO splittings AE-p^o can be directly un- 
derstood by the PSS breaking term ly^^^ within the representation of the SUSY partner 
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Hamiltonian H. Furthermore, this symmetry breaking term can be treated as a very small 
perturbation on the exact PSS limit Hq^^- Therefore, the perturbative nature of PSS can 
be pinned down. 

IV. SUMMARY AND PERSPECTIVES 

The supersymmetric (SUSY) quantum mechanics is used to investigate the origin of the 
pseudospin symmetry (PSS) and its breaking in the single-particle spectra in nuclei. In the 
formalism, it can be seen that, while the spin symmetry conserving term k,{k + 1) appears 
in the single-particle Hamiltonian H, the PSS conserving term k{k, — 1) appears naturally 
in its SUSY partner Hamiltonian H. In addition, the fact that all states with / > have 
their own pseudospin partners but not the so-called intruder states can be interpreted by 
employing both exact and broken SUSY in a unified way. 

In the present study, we focus on a Schrodinger equation with a Woods-Saxon central 
potential. This corresponds to the lowest-order approximation for transforming a Dirac 
equation into a diagonal form by using the similarity renormalization group (SRG) technique. 

Within the single-particle Hamiltonian H representation, the wave functions of pseu- 
dospin doublets are very different from each other, and the small pseudospin-orbit (PSO) 
splittings AEpso are due to the sophisticated cancellations among all channels. However, 
this does not necessarily mean the dynamical or even the nonperturbative nature of PSS. 
By taking the Hamiltonian with a harmonic oscillator (HO) potential H^'~' as the exact PSS 
limit, it is found that the largest perturbation correction due to the symmetry breaking term 
y^HO jg roughly 0.13, and Ai^pso can be well reproduced by the third-order perturbation 
calculations. These results pin down the perturbative nature of PSS. 

The origin of the PSS and its breaking mechanism can be interpreted explicitly within the 
representation of the SUSY partner Hamiltonian H. In order to obtain H, one first solves 
the first-order differential equation for the reduced superpotentials ^'^(r), and then obtains 
the K-dependent central potentials V^{r) in H. It is important that the central potentials 
are regular and converge at both r — and r — oo. Consequently, the eigenstates of 
Hamiltonians H and H are exactly one-to-one identical, except the additional eigenstates 
in H when the SUSY is exact, which correspond to the intruder states without pseudospin 
partners. 
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By holding this one-to-one mapping relation, the origin of PSS deeply hidden in H can be 
traced by employing its SUSY partner Hamiltonian H. Within this SUSY partner scheme, 
the wave functions of the pseudospin doublets are almost identical, and the well known node 
relation between the pseudospin doublets, i.e., n = n — 1 for n < and n = n for k > 0, 
is one of the intrinsic properties of the SUSY quantum mechanics. It is also found that the 
PSO splittings AE'pso are mainly contributed by the PSS breaking terms. This indicates 
the sophisticated cancellations among different terms in H can be understood in a clear 
scheme by using a proper decomposition. 

The Hamiltonian H is then split into the PSS conserving term Hq^^ and the PSS breaking 
term 14^^^^. While the PSS conserving potentials Vpss('") remain an approximate Woods- 
Saxon shape with a slight /t-dependence, the PSS breaking potentials V^Q,o{r) show several 
special features. 1) The potentials Vpa,o{r) are regular functions of r, and vanish at both 
r — and r — t- oo. 2) The amplitudes of Vpso are around 1 MeV, which directly lead to 
the reduced PSO splittings A^Epso ^ 1 MeV. 3) The potentials Vpso('^) show a particular 
shape that they are negative at small radius but positive at large radius, and with a node 
at the surface region. The general pattern that the PSO splittings become smaller with the 
increasing single-particle energies, and even reverse for resonance states, can be understood 
straightforwardly by such particular symmetry-breaking potentials. 

Finally, the perturbation calculations are performed based on the PSS conserving Hamil- 
tonian Hq^^ with the perturbation VT^^^. It is found that in such decomposition the largest 
perturbation correction due to the symmetry breaking term VT^^^ is less than 0.03, and not 
only the PSO splittings but also the energy degeneracy of the spin doublets are excellently 
reproduced by the first-order perturbation calculations. In such way, the origin of the PSS 
can be recognized, and its breaking is due to a very small perturbation 14^^^^ on the exact 
PSS limit ^o^ss 

The present investigation employs a Schrodinger equation for illustrating the key ideas on 
applying the SUSY quantum mechanics to the PSS in nuclei. As one of the natural exten- 
sions, a remarkable spin-orbit potential should also be implemented for a more appropriate 
description of the single-particle states. In order to completely answer the question why 
the PSS is conversed better than the SS in realistic nuclei, the intrinsic relation between 
the spin-orbit potential and the central potential or the effective mass must be taken into 
account. In this sense, the PSS must be regarded as the relativistic symmetry, and it should 
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be recognized in the Dirac equation, or equivalently the Schrodinger-hke equation obtained 
by transforming the Dirac equation into a diagonal form by using the SRG technique. 
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